Genomic characterization of human respiratory syncytial virus circulating in Islamabad, Pakistan, during an outbreak in 2022-2023.

In this study, conducted at the National Institute of Health, Islamabad, during an outbreak of human respiratory syncytial virus (hRSV) from December 2022 to January 2023, the first whole-genome sequences of hRSV isolates from Islamabad, Pakistan, were determined. Out of 10 positive samples, five were sequenced, revealing the presence of two genotypes: RSV-A (GA2.3.5, ON1 strain) and RSV-B (GB5.0.5.a, BA-10 strain). A rare non-synonymous substitution (E232G) in G the protein and N276S in the F protein were found in RSV-A. In RSV-B, the unique mutations K191R, Q209R, and I206M were found in the F protein. These mutations could potentially influence vaccine efficacy and viral pathogenicity. This research underscores the importance of genomic surveillance for understanding RSV diversity and guiding public health responses in Pakistan.

Insights into measles virus: Serological surveillance and molecular characterization.

BACKGROUND: Measles has been a significant public health concern in Pakistan, especially in the Khyber Pakhtunkhwa (KPK) province, where sporadic and silent epidemics continue to challenge existing control measures. This study aimed to estimate the prevalence and investigate the molecular epidemiology of the measles virus (MeV) in KPK and explore the vaccination status among the suspected individuals. METHODS: A cross-sectional study was conducted between February and October 2021. A total of 336 suspected measles cases from the study population were analyzed for IgM antibodies using Enzyme-Linked Immunosorbent Assay (ELISA). Throat swabs were randomly collected from a subset of positive cases for molecular analysis. Phylogenetic analysis of MeV isolates was performed using the neighbor-joining method. The vaccination status of individuals was also recorded. RESULTS: Among the suspected participants, 61.0% (205/336) were ELISA positive for IgM antibodies, with a higher prevalence in males (64.17%) compared to females (57.04%). The majority of cases (36.0%) were observed in infants and toddlers, consistent with previous reports. The majority of IgM-positive cases (71.7%) had not received any dose of measles vaccine, highlighting gaps in vaccine coverage and the need for improved immunization programs. Genetic analysis revealed that all MeV isolates belonged to the B3 genotype, with minor genetic variations from previously reported variants in the region. CONCLUSION: This study provides valuable insights into the genetic epidemiology of the MeV in KPK, Pakistan. The high incidence of measles infection among unvaccinated individuals highlights the urgency of raising awareness about vaccine importance and strengthening routine immunization programs.

Crimean-Congo Hemorrhagic Fever Virus Diversity and Reassortment, Pakistan, 2017-2020.

Sporadic cases and outbreaks of Crimean-Congo hemorrhagic fever (CCHF) have been documented across Pakistan since 1976; however, data regarding the diversity of CCHF virus (CCHFV) in Pakistan is sparse. We whole-genome sequenced 36 CCHFV samples collected from persons infected in Pakistan during 2017-2020. Most CCHF cases were from Rawalpindi (n = 10), followed by Peshawar (n = 7) and Islamabad (n = 4). Phylogenetic analysis revealed the Asia-1 genotype was dominant, but 4 reassorted strains were identified. Strains with reassorted medium gene segments clustered with Asia-2 (n = 2) and Africa-2 (n = 1) genotypes; small segment reassortments clustered with the Asia-2 genotype (n = 2). Reassorted viruses showed close identity with isolates from India, Iran, and Tajikistan, suggesting potential crossborder movement of CCHFV. Improved and continuous human, tick, and animal surveillance is needed to define the diversity of circulating CCHFV strains in Pakistan and prevent transmission.

Evolutionary analysis of seasonal influenza A viruses in Pakistan 2020-2023.

INTRODUCTION: Influenza A viruses cause global health concerns due to their high amino acid substitution rates. They are linked to yearly seasonal epidemics and occasional pandemics. This study focused on sequencing influenza A virus strains in Pakistan. MATERIALS AND METHODS: We analyzed the genetic characteristics of influenza A(H1N1)pdm09 and A(H3N2) viruses circulating in Pakistan from January 2020 to January 2023. Whole genome sequences from influenza A (n = 126) virus isolates were amplified and sequenced by the Oxford Nanopore (MinION) platform. RESULTS: The HA genes of influenza A(H1N1)pdm09 underwent amino acid substitutions at positions K54Q, A186T, Q189E, E224A, R259K, and K308R in sequenced samples. The HA genes of influenza A(H3N2) had amino acid substitutions at G53D, E83K, D104G, I140M, S205F, A212T, and K276R in the sequenced samples. Furthermore, the HA gene sequences of influenza A(H1N1)pdm09 in this study belonged to subclade 6B.1A.5a.2a. Similarly, the HA gene sequences of influenza A(H3N2) were classified under six subclades (3C.3a.1 and 3C.2a1b.2a [2, 2a.1, 2b, 2c, and 2a.3b]). Notably, amino acid substitutions in other gene segments of influenza A(H1N1)pdm09 and A(H3N2) were also found. CONCLUSION: These findings indicate influenza A(H1N1)pdm09 and A(H3N2) viruses co-circulated during the 2020-2023 influenza season in Pakistan. Continued surveillance is crucial for real-time monitoring of possible high-virulence variation and their relevance to existing vaccine strains.

Integrated approach for detection of SARS-CoV-2 and its variant by utilizing LAMP and ARMS-PCR.

Global impact of COVID-19 pandemic has heightened the urgency for efficient virus detection and identification of variants such as the Q57H mutation. Early and efficient detection of SARS-CoV-2 among densely populated developing countries is paramount objective. Although RT-PCR assays offer accuracy, however, dependence on expansive kits and availability of allied health resources pose an immense challenge for developing countries. In the current study, RT-LAMP based detection of SARS-Cov-2 with subsequent confirmation of Q57H variant through ARMS-PCR was performed. Among the 212 collected samples, 134 yielded positive results, while 78 tested negative using RT-LAMP. Oropharyngeal swabs of suspected individuals were collected and processed for viral RNA isolation. Isolated viral RNA was processed further by using either commercially available WarmStart Master Mix or our in house developed LAMP master mix separately. Subsequently, the end results of each specimen were evaluated by colorimetry. For LAMP assays, primers targeting three genes (ORF1ab, N and S) were designed using PrimerExplorer software. Interestingly, pooling of these three genes in single reaction tube increased sensitivity (95.5%) and specificity (93.5%) of LAMP assay. SARS-CoV-2 positive specimens were screened further for Q57H mutation using ARMS-PCR. Based on amplicon size variation, later confirmed by sequencing, our data showed 18.5% samples positive for Q57H mutation. Hence, these findings strongly advocate use of RT-LAMP-based assay for SARS-CoV-2 screening within suspected general population. Furthermore, ARMS-PCR also provides an efficient mean to detect prevalent mutations against SARS-Cov-2.

Genomic characterization of the first mpox and varicella-zoster co-infection in Pakistan (2023) through next-generation sequencing.

The global mpox outbreak spanning 2022-2023 has affected numerous countries worldwide. In this study, we present the first report on the detection, whole-genome sequence, and coinfection of the mpox virus and varicella zoster virus (VZV) from Pakistan. During April-May 2023, samples from 20 suspected cases of mpox were tested at the National Institutes of Health, Islamabad among which 4 tested positive. All four cases had a travel history of Saudi Arabia. All the suspected samples were processed by using a Zymo research kit for DNA extraction, followed by qRT-PCR amplification by using a DaAn Gene detection kit for the mpox virus. Further, two of the positive samples with a low Ct value (<20) were subjected to whole-genome sequencing using a metagenomic approach on the iSeq (Illumina) platform. The sequencing results revealed Clade IIb and genotype A.2.1 of MPXV, which clustered with viruses from Slovenia and the UK in July and June 2022, respectively. Our analysis identified two novel nonsynonymous substitutions in mpox virus, namely V98I in OPG046 and P600S in OPG109. Furthermore, we successfully retrieved the complete genome of VZV from the same sample, belonging to Clade 5. This study represents the first positive case of MPXV in Pakistan and the coinfection of mpox and VZV by using a metagenome approach providing insights into their complete genomes. Our results highlight the importance of surveillance at the point of entries, strengthening lab capacities including next-generation sequencing, and using differential diagnosis for timely and accurate detection of mpox cases.

Whole-genome sequencing of Crimean-Congo hemorrhagic fever virus circulating in Pakistan during 2022.

Pakistan is an endemic country for Crimean-Congo hemorrhagic fever (CCHF) and its Balochistan province is considered a hotspot region for circulation of the virus whereas sporadic cases have been reported from other parts of the country. Our study aims to investigate the genomic diversity of the CCHF virus circulating in Punjab and Khyber Pakhtunkhwa provinces of Pakistan. Between April to September 2022, 46 samples from suspected CCHF patients were tested, with 6 (13%) showing positive RT-PCR results. Among the positive cases, all were male, aged 14-48 years among which 4 were butchers. Three CCHF patients succumbed to the disease. The complete S-M-L-gene fragments of 4 positive samples were sequenced. The S and L segments belonged to the Asia-1 clade and clustered with regional strains from Iran, India, and Afghanistan. One M segment sequence grouped into Africa-2 along with those reported from India during 2016-2019. We report the detection of a reassorted virus (NIH-PAK-CCHF-03|2022) having Asia-1-Africa-2-Asia-1 (S-M-L) segment pattern for the first time from Pakistan. Mutational analysis showed M segments harboring eight mutations (T55A, S80P, T110I, T185A, T189A, A212T, and N239I/T) in the mucin-like domain, five mutations (D250N, T333S, I375V, M401I, A433T), four mutations (N545D, Y657F, K688R, and I824V) in GP38-domain, and three mutations (T1418N, A1431V, and G1449S) in Gc-domain. These findings highlight the significance of whole-genome sequencing of indigenous strains for a better understanding of the CCHFV evolution in Pakistan.

First report of coinfection and whole-genome sequencing of norovirus and sapovirus in an acute gastroenteritis patient from Pakistan.

Acute gastroenteritis is one of the most common diseases in infants and children in developing countries including Pakistan. In Pakistan, rotavirus (RVA) is known to contribute significantly to pediatric diarrheal illness, but the contribution of other viruses is still unclear. In the current study we have identified a case of mixed infection of norovirus (NoV) and sapovirus (SaV) in a 2-year-old child with acute gastroenteritis. The sample was initially processed for the detection of group A RVA through ELISA followed by NoV using RT-PCR assay. The sample tested positive for NoV RNA and was later subjected to whole-genome sequencing using meta-genome approach on Miseq (Illumina) platform. Sequencing results revealed GII.15 genotype of NoV that clustered with viruses from China and USA from 2017 to 2021. We also retrieved the complete genome of SaV (GI.1 genotype) from the same sample and phylogenetic analysis showed clustering with strains reported from Japan, South Korea, US, and Taiwan during 2012-2016. This is the first report from Pakistan that confirms coinfection of NoV and SaV and elucidates their whole genomes. We recommend initiation of NoV and SaV surveillance program to ascertain disease burden and explore genetic diversity, especially as RVA vaccines have been included in national immunization program.

Genomic Characterization of Dengue Virus Outbreak in 2022 from Pakistan.

Pakistan, a dengue-endemic country, has encountered several outbreaks during the past decade. The current study aimed to explore the serotype and genomic diversity of dengue virus responsible for the 2022 outbreak in Pakistan. From August to October 2022, NS-1 positive blood samples (n = 343) were collected from dengue patients, among which, (85%; n = 293) were positive based on RT-PCR. In terms of gender and age, dengue infection was more prevalent in male patients (63%; n = 184), with more adults (21-30 years; n = 94) being infected. The serotyping results revealed DENV-2 to be the most predominant serotype (62%; n = 183), followed by DENV-1 (37%; n = 109) and DENV-3 (0.32%; n = 1). Moreover, a total of 10 samples (DENV-2; n = 8, DENV-1; n = 2) were subjected to whole-genome sequencing. Among these, four were collected in early 2022, and six were collected between August and October 2022. Phylogenetic analysis of DENV-2 sequenced samples (n = 8) revealed a monophyletic clade of cosmopolitan genotype IVA, which is closely related to sequences from China and Singapore 2018, and DENV-1 samples (n = 2) show genotype III, which is closely related to Pakistan isolates from 2019. We also reported the first whole genome sequence of a coinfection case (DENV1-DENV2) in Pakistan detected through a meta-genome approach. Thus, dengue virus dynamics reported in the current study warrant large-scale genomic surveillance to better respond to future outbreaks.

Breakthrough cases of Omicron and Delta variants of SARS-CoV-2 during the fifth wave in Pakistan.

COVID-19 pandemic has severely affected Pakistan with 1,557,134 cases as of August 4, 2022. However, the data regarding breakthrough infections in Pakistan is scant. Hence, the objective was to analyze SARS-CoV-2 breakthrough infections with respect to vaccines and variants during the fifth wave in Pakistan. Therefore, the Department of Virology (NIH, Pakistan) genotyped 2,467 randomly selected individuals between November 2021 and February 2022 using the SNPsig® SARS-CoV-2 (EscapePLEX) kit (PrimerDesign, UK). P681R and K417N mutations were used to distinguish delta and omicron. Data on the patient's age, gender, date of collection, variant, and vaccination status were analyzed using Statistical Package for Social Sciences (SPSS) software. Among 2,467 genotyped samples, Omicron was detected in 58.6% (n = 1445), Delta in 40.4% (n = 998) and undetermined/wildtype variant in 24 samples. The vaccination status of omicron-positive patients showed (49.7%; n = 718/1445) and Delta-positive patients (39.67%; n = 396/998) to be fully vaccinated. Of note, a high percentage 85% of breakthrough cases (n = 947) were identified among fully vaccinated individuals (n = 1114). Among them, 85.9% (n = 617/718) belonged to omicron and 83.3% (n = 330/396) to delta. Moreover, 76.7% (n = 855) of vaccinated individuals (n = 1114) received Sinopharm (n = 432) and Sinovac (n = 423) vaccines. The majority of breakthrough subjects who contracted Omicron were vaccinated with Sinopharm (93.0%, n = 256) and delta with Cansino (100%, n = 44). Individuals vaccinated with Sinovac showed the most frequent breakthrough cases for both Omicron and Delta variant between the 4th and 6th months (n = 278) after primary vaccination as compared to the 7th to 9th months (n = 24) category. While in case of Sinopharm, maximum breakthrough cases occurred between 7th to 9th months (n = 234) as compared to the 4th to 6th months (n = 120) after primary vaccination. Omicron and Delta breakthrough cases in men (n = 364 and 193) are more frequently seen than women (n = 253 and 138) respectively and breakthrough majority cases (n = 392) occurred in individuals aged 18-33 years. Breakthrough cases limiting monitoring in Pakistan impose a substantial constraint on policymakers' ability to take timely effective decisions. Since the current study consists of only a 2,467-genotyped sample, comprehensive data should be analyzed.

Genomic diversity of SARS-CoV-2 in Pakistan during the fourth wave of pandemic.

The emergence of different variants of concern of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in upsurges of coronavirus disease 2019 (COVID-19) cases around the globe. Pakistan faced the fourth wave of COVID-19 from July to August 2021 with 314,786 cases. To understand the genomic diversity of circulating SARS-CoV-2 strains during the fourth wave of the pandemic in Pakistan, this study was conducted. The samples from 140 COVID-19-positive patients were subjected to whole-genome sequencing using the iSeq Sequencer by Illumina. The results showed that 97% (n = 136) of isolates belonged to the delta variant while three isolates belonged to alpha and only one isolate belonged to the beta variant. Among delta variant cases, 20.5% (n = 28) isolates were showing B.1.617.2 while 23.5% (n = 25), 17.59% (n = 19), 14.81% (n = 16), and 13.89% (n = 15) of isolates were showing AY.108, AY.43 AY.127, and AY.125 lineages, respectively. Islamabad was found to be the most affected city with 65% (n = 89) of delta variant cases, followed by Karachi (17%, n = 23), and Rawalpindi (10%, n = 14). Apart from the characteristic spike mutations (T19R, L452R, T478K, P681R, and D950N) of the delta variant, the sublineages exhibited other spike mutations as E156del, G142D, T95I, A222V, G446V, K529N, N532S, Q613H, and V483A. The phylogenetic analysis revealed the introductions from Singapore, the United Kingdom, and Germany. This study highlights the circulation of delta variants (B.1.617.2 and sublineages) during the fourth wave of pandemic in Pakistan.

Molecular identification and characterization of measles virus genotypes circulating in Afghanistan, 2008-2018.

The measles virus (MV) remains a leading cause of morbidity and mortality in children under 5 years of age. Molecular identification of circulating wild-type MV) strains is a vital component of the measles elimination program. We received 159 oral swab samples from Afghanistan during 2008-2018. Viral RNA was extracted, followed by one-step RT-PCR and positive amplicons were subject to sequencing for genotype identification. Out of 159 total samples, 52% (83/159) were detected positive by RT-PCR. Genotype D4 was identified from 2.4% (2/83), genotype H1, 4.8% (4/83), and genotype B3, 92.7% (77/83) cases, respectively.

Molecular detection and genetic characterization of human metapneumovirus strains circulating in Islamabad, Pakistan.

Lower respiratory illness is one of the leading causes of death among children in low- and high-income countries. Human metapneumovirus (hMPV) is a key contributor to respiratory illnesses commonly reported among children and causes serious clinical complications ranging from mild respiratory infections to severe lower respiratory tract anomalies mainly in the form of bronchiolitis and pneumonia. However, due to the lack of a national surveillance system, the clinical significance of hMPV remains obscure in the Pakistani population. This study was conducted to screen throat swabs samples collected from 127 children reported with respiratory symptoms at a tertiary care hospital in Islamabad. Out of 127, 21 (16.5%) samples were positive for hMPV with its genotype distribution as A2a (10%), A2b (20%), B1 (10%), and B2 (60%). Phylogenetic analysis showed that the hMPV viruses were closely related to those reported from neighboring countries including India and China. This work will contribute to a better understanding of this virus, its diagnosis, and the handling of patients in clinical setups. Further studies at a large-scale are warranted for a better understanding of the disease burden and epidemiology of hMPV in Pakistan.